The role of local and global strangeness neutrality at the inhomogeneous freeze-out in 

relativist ic heavy ion collisions 
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The decoupling surface in relativistic heavy-ion collisions may not be homogeneous. Rather, 
inhomogeneities should form when a rapid transition from high to low entropy density occurs. 
We analyze the hadron "chemistry" from high-energy heavy-ion reactions for the presence of such 
density inhomogeneities. We show that due to the non-linear dependence of the particle densities on 
the temperature and baryon-chemical potential such inhomogeneities should be visible even in the 
integrated, inclusive abundances. We analyze experimental data from Pb+Pb collisions at CERN- 
SPS and Au+Au collisions at BNL-RHIC to determine the amplitude of inhomogeneities and the 
role of local and global strangeness neutrality. 



INTRODUCTION 

It is expected that at sufficiently high energies, a 
transient state of deconfined matter with broken Z(3) 
center symmetry and/or with (approximately) restored 
chiral symmetry is produced in collisions of heavy nu- 
clei. Lattice QCD simulations [2] indicate that a second- 
order critical point exists, which was predicted by effec- 
tive chiral Lagrangians [3]; present estimates locate it 
at T ss 160 MeV, \x B 360 MeV. This point, where 
the cr-ficld is massless, is commonly assumed to be the 
cndpoint of a line of first-order phase transitions in the 
(fiB,T) plane. To detect that endpoint, it is hoped 
that by varying the beam energy, for example, one can 
"switch" between the regimes of first-order phase transi- 
tion and cross over, respectively. If the particles decouple 
shortly after the expansion trajectory crosses the line of 
first order transitions one may expect a rather inhomo- 
geneous (energy-) density distribution on the freeze-out 
surface [4, 5] (similar, say, to the CMB photon decoupling 
surface observed by WMAP [6]). On the other hand, if 
the low-temperature and high-temperature regimes are 
smoothly connected, pressure gradients tend to wash out 
density inhomogeneities. Similarly, in the absence of 
phase-transition induced non-equilibrium effects, the pre- 
dicted initial-state density inhomogeneities [7, 8] should 
be strongly damped. 

Thus, we investigate the properties of an inhomoge- 
neous fireball at (chemical) decoupling. Note that if the 
scale of these inhomogeneities is much smaller than the 
decoupling volume then they can not be resolved indi- 
vidually, nor will they give rise to large event-by-event 
fluctuations. Because of the nonlinear dependence of the 
hadron densities on T and [is, they should nevertheless 
reflect in the average abundances. 

Our basic assumption is that as the fireball expands 
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and cools, at some stage the abundances of hadrons 
"freeze" , keeping memory of the last instant of chemical 
equilibrium. This stage is refered to as chemical freeze- 
out. By definition, only processes that conserve particle 
number for each species individually, or decays of unsta- 
ble particles may occur later on. The simplest model is 
to treat the gas of hadrons within the grand canonical 
ensemble, assuming a homogeneous decoupling volume. 
The abundances are then determined by two parameters, 
the temperature T and the baryonic chemical potential 
fiB] the chemical potential for strangeness is fixed by 
the condition for overall strangeness neutrality. Fits of 
hadronic ratios were performed extensively [9, 10] within 
this model, sometimes also including a strangeness (7 S ) 
or light quark (j q ) supression factor [11, 12] or interac- 
tions with the chiral condensate [13]. 

In [14] we analyzed the experimental data on relative 
abundances of hadrons with respect to the presence of in- 
homogeneities on the decoupling surface. To that end we 
proposed a very simple and rather schematic extension of 
the common grand canonical freeze-out model, i.e. a su- 
perposition of such ensembles with different temperatures 
and baryon-chemical potentials. Each ensemble is sup- 
posed to describe the local freeze-out on the scale of the 
correlation length ~ 1/T ~ 1 — 2 fm. Even if freeze-out 
occurs near the critical point, the correlation length of the 
chiral condensate is bound from above by finite size and 
finite time effects, effectively resulting in similar num- 
bers [15]. On the other hand, for small chemical poten- 
tial, far from the region where the cr-field is critical, the 
relevant scale might be set by the correlation length for 
Polyakov loops, which is of comparable magnitude [16]. 
Classical nucleation theory for strong first-order phase 
transitions predicts even larger "bubbles" [17] but is un- 
likely to apply to small, rapidly expanding systems en- 
countered in heavy-ion collisions [5, 18]. Another (classi- 
cal) model for the formation of small droplets in rapidly 
expanding QCD matter has been introduced in [19]. The 
entire decoupling surface contains many such "domains" , 
even if a cut on mid-rapidity is performed. We therefore 
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expect that the distributions of temperature and chem- 
ical potential are approximately Gaussian [20]. Besides 
simplicity, another goal of the present analysis is to avoid 
reference to a particular dynamical model for the forma- 
tion or for the distribution of density perturbations. In 
fact, we presently aim merely at checking whether any 
statistically significant signal for the presence of inhomo- 
geneities is found in the data. If so, more sophisticated 
dynamical models could be employed in the future to 
understand the evolution of inhomogeneitics from their 
possible formation in a phase transition until decoupling. 

Rate equations for nuclear fusion and dissociation pro- 
cesses (and neutron diffusion) have been used for inho- 
mogeneous big bang nucleosynthesis in the early uni- 
verse [21]. Similarly, hadronic cascade models could be 
used for heavy-ion reactions [22] . This would remove ref- 
erence to the grand-canonical ensemble and to a thin de- 
coupling surface in space-time. In fact, hadronic binary 
rescattering models do predict a rather thick freeze-out 
layer [22, 23], where matter expands non- ideally. On the 
other hand, the steep drop of multi-particle collision rates 
with temperature should narrow the freeze-out again [24] . 
In either case, we do not expect a strong energy depen- 
dence of the width of freeze-out (see also [25] ) . 

At chemical freeze-out, matter is in a state of expan- 
sion. However, such flow effects do not affect the rela- 
tive abundances of the particles (in full phase space) if 
their densities are homogeneous throughout the decou- 
pling volume. The total number of particles of species i, 
integrated over a solid angle of 4tt, is given by an integral 
of the current Nf* = piu^ , with the four- velocity of 
the expanding fluid, over a given freeze-out hypersurface 

o» = (t fo ,x io y. 



N, t = I da,N? = Pl (T l °,^) 



(1) 



The second factor on the r.h.s. is nothing but the three- 
volume V3 of the decoupling hypersurface as seen by the 
observer. This volume is common to all species and drops 
out of multiplicity ratios: Ni/Nj = p f °/p f °. It is clear 
that the argument holds even when cuts in momentum 
space are performed, provided that the differential distri- 
butions of all particles do not depend on that particular 
momentum-space variable (for example, rapidity cuts for 
boostinvariant expansion [26]). 

When the intensive variables T and [is vary, then the 
integration measure (/ u ■ da) /V3 will, in general, depend 
on the assumed distribution and amplitude of inhomo- 
geneities, as well as on the hydrodynamic flow profile etc. 
Nevertheless, it is still the same for all particle species 
and so can be written in the form 



Ju-da — ► J dTdp B P(T, p B ) 



(2) 



with P(T,/j,b) some distribution for T and /ig. For sim- 
plicity, and for lack of an obvious motivation for assum- 
ing otherwise, we shall take P(T,/j,b) to factorize into 



a distribution for T, times one for /is- These distribu- 
tions could, in principle, be obtained from the real-time 
evolution of the phase transition [4, 5]. 

THE MODEL 

In [14] we introduced our model to analyze the avail- 
able data from heavy-ion collisions at CERN-SPS and 
BNL-RHIC. There the hadron abundances are deter- 
mined by four parameters: the arithmetic means of the 
temperatures and chemical potentials of all domains, T 
and ~p Bl and the widths of their Gaussian distributions, 
ST and Sfi B - Of course, the densities of strange par- 
ticles depend also on the strangeness-chemical poten- 
tial [is, which we determined in [14] by imposing local 
strangeness neutrality. That means, the strange chem- 
ical potential in each single domain was fixed by de- 
manding zero net strangeness there. However, the effect 
of independent fluctuations of ps should also be looked 
at, in particular for collisions at low and intermediate 
energies (^/snn^ 15 GcV). This may help for example 
to reproduce the A to p ratio, which was found to be 
larger than one [28] and the K + /tt + enhancement around 
E~Lat>/A — 30 GeV [29]. Allowing for such independent 
fluctuations, the hadron abundances depend on six pa- 
rameters: the arithmetic means of the temperatures and 
chemical potentials of all domains, T,~fi B and Ji s , and 
the widths of their Gaussian distributions, ST,S[ib and 
5 [is- They read: 

oc 

~ Pl (T,JI B ,JI s ,6T,6[iB,6[is) = J dTP(T;T,ST) (3) 



00 00 

J d[i B P{[ib;~Pb,Svb) j d[i s P(fi,s;~Ps> s l J 's) Pi(T,[i B ,[is) , 

— 00 —00 

with pi(T, [iB, Ms) the actual "local" density of species i, 
and with P(x;x,Sx) <~ cxp[— (x — x) 2 /2 5x 2 ] the distri- 
bution of temperatures and chemical potentials within 
the decoupling three-volume (the proportionality con- 
stants normalize the distributions over the intervals 
where they are defined). In addition, strangeness con- 
servation enters now as a global constraint for the mean 
of the strange chemical potential JI S : 

fs = YjPi (T^B^s^ T ^^B,S[is)(nl~n^) = 0, (4) 

i 

with f s the net-strangeness, n\ , rij the number of strange 
and anti-strange quarks of hadron species i, respectively. 
That means, the global densities obtained for given val- 
ues of temperature and chemical potential parameters 
weighted with the corresponding net number of strange 
quarks are summed and demanded to vanish to guaran- 
tee strangeness neutrality. In the limit ST, Sps, Sps —> 
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the Gaussian distributions are replaced by <5-functions 
and the conventional homogeneous freeze-out scenario is 
recovered: 

Pi (T,7* B , 7* s , 0,0,0) = pi(T,Jl B ,-p s ), (5) 

and the corresponding strangeness neutrality condition 
fixing ~pZ s . In other words, in that limit the average den- 
sities are uniquely determined by the first moments of 
T and ps- For the present investigation, we set the 
width of the distribution for the strange chemical po- 
tential equal to zero, Sps = 0. Since eq.(4) only ensures 
global strangeness neutrality, in this limit still finite net 
strangeness values in individual domains will appear, in 
contrast to our former analysis, where we fixed p s by 
f s = locally. It is important to note that with setting 
Sp s = and the global constraint eq. 4 for ~ps, the den- 
sities again are a function of four parameters: T,~jI B ,ST 
and S/ib- Thus we will write all quantities again as a 
function of these four parameters only. In the follow- 
ing we will investigate how the fits to the experimentally 
measured particle abundances arc influenced by the dif- 
ferent strangeness neutraliy conditions [31]. 

For the present analysis we compute the densities 
Pi{T,ps) in the ideal gas approximation, supplemented 
by an "excluded volume" correction: 

MT,M = TT -^-^. (6) 

This schematic correction models repulsive interactions 
among the hadrons at high densities. Vi denotes the 
volume occupied by a hadron of species i; we employ 
v = f7r.R0 3 with R = 0.3 fm for all species [27]. There- 
fore, for the homogeneous model the denominator in (6) 
drops out of multiplicity ratios. This is not the case for 
an inhomogeneous decoupling surface, where the distri- 
butions of various species differ. For all fits over the full 
solid angle, we fixed the isospin chemical potential by 
equating the total charge in the initial and final states; 
for the mid-rapidity fits at high energies, we fixed pi = 0. 

To illustrate the effect of inhomogeneities on the distri- 
butions of various hadrons within the decoupling volume 
we introduce 

A (T; T, p B , ST, Sp B )=P(T; T, ST) 

OO 

J dp B P(p B ;p B ,Sp B ) p l (T,p B ) 

x- = , (7) 

pAT^b.ST^pb) 

Di(p B ;T,JI B ,ST,6p B ) = P(p B ;p B ,Sp B ) 

00 

/ dT P(T; T, ST) Pi (T,p B ) 

x ° _ __ . (8) 

ft (T,p B -,ST,Sp B ) 

Di(T), for example, is the probability that a particle of 
type i was emitted from a domain of temperature T. 



The main contribution to the integrals in (3) is not from 
T and ~p B since hot spots shine brighter than "voids". 
Rather, they are dominated by the stationary points of 
the distributions defined in eqs. (7,8) above. Hence, the 
average emission temperature (T% and baryon-chemical 
potential {pb)i in general depend on the particle species 
i, unless ST = Sps = 0. They can be evaluated as 

oc 

(T)i = J dTT Di(T;T,JiB,6T,S»B) , 


00 

(p B )i = J dp B Pb Di(pB;T 7 Ji B7 ST,SpB) ■ (9) 

—00 

Physically, this means that for non-zero widths of the 
temperature and chemical potential distributions the 
freeze-out volume is not perfectly "stirred", in that the 
relative concentrations of the particles vary. 

DATA ANALYSIS 

To determine the four parameters of the model we min- 
imize 

x 2 = J2^ p -rr del ) 2 /^ (10) 

i 

in the space of T, 7i B , ST, and Sps- That is, we obtain 
least-square estimates for the parameters, assuming that 
they are independent. In (10), r^ xp and r" Lodel denote 
the experimentally measured and the calculated particle 
ratios, respectively, and of is set by the uncertainty of the 
measurement. Wherever available, we sum systematic 
and statistical errors in quadrature. 

The data used in our analysis are the particle multi- 
plicities measured by the NA49 collaboration for central 
Pb+Pb collisions at beam energy E Lah /A = 20, 30, 40, 80 
and 158 GeV [29], and those measured by STAR for cen- 
tral Au+Au collisions at BNL-RHIC, ref. [30] {^/snn = 
130 GeV, compiled in [32]) and ref. [33] (200 GeV). At 
RHIC energies, we analyze the midrapidity data; at top 
SPS energy, both, midrapidity and 47r data. At all other 
energies, we restrict ourselves to the An solid angle data 
by NA49 in order to avoid biases arising from differing 
acceptance windows of various experiments. Further- 
more, our checks showed that the fit results can depend 
somewhat on the actual selection of experimental ratios. 
Hence, where possible, we have opted for the least bias 
by choosing r^ xp — N^ xp /N% xp , i.e. the multiplicity of 
species i relative to that of pions. This represents the 
maximal set of independent data points, as it is equiv- 
alent to fitting absolute multiplicities with an additional 
overall three- volume parameter, Ni = Vspi. 

Specifically, at E Lah /A = 20, 30, and 80 GeV the mul- 
tiplicities of 7T+, 7r~, K + , K~, B — B, A, A, and <f) are 
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available. For the (in-)homogeneous model, this leaves 
five (three) degrees of freedom. At 40 GeV, we can add 
the H~ and fi + fi. The data sets for top SPS energies in- 
clude yet a few more species: p, p (only midrapidity), Kg 
(only 47r), S + and fl, fl seperately. For RHIC-130, we 
fitted to the K+/K~, p/p, A/A, 3+/S - , U/Q, K'/tt', 
K s /ir~, p/tt~, A/7r~, Kq/tt~, <\>/-k~ , E~ /tt~ and f2/7r~ 
ratios. Finally, at RHIC-200 the K+/K~, p/p, Q/fi, 

K~ /lT~, p/Tr~, A/tT~, A/7T~, S^/tT^, S + /7T~, O / ir~ , 

4>/K~ and Kq/K~ ratios were used. The first three ra- 
tios are close to unity and essentially just set the chemical 
potentials to zero; they do not help to fix T, ST and SfiB- 
Where appropriate, feeding from strong and electro- 
magnetic decays has been included in r ™' odel by replacing 
Pi — > Pi + Bij pj . The implicit sum over j i runs over 
all unstable hadron species, with Bij the branching ratio 
for the decay j — ► i, which were taken from [34]. From 
all the resonances listed by the Particle Data Group [34] , 
mesons up to a mass of 1.5 GeV and baryons up to a mass 
of 2 GeV were included, respectively. The finite widths 
of the resonances were not taken into account, and un- 
known branching ratios were excluded from the feeding. 
These details are irrelevant for the qualitative behavior 
of ST and Sps but do, of course, matter for quantitative 
results. 

RESULTS 
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FIG. 1: X 2 /dof versus t/snn for the homogeneous (ST = 
S/j, — 0, squares) and the inhomogeneous fit (ST and S/j, free 
parameters, circles and triangles). Circles denote the case 
of local strangeness neutrality, while triangles represent the 
global strangeness neutrality case. The lines are meant to 
guide the eye. Furthermore, the x 2 /dof corresponding to the 
95.4% confidence interval is shown by the dotted line. 

Fig. 1 shows the minimal x 2 P er degree of freedom 
(taken as the number of data points minus the number 



of parameters) for the homogeneous approach and the in- 
homogeneous approach with local or global strangeness 
neutrality, respectively. Note that the x 2 -values for the 
homogeneous model are in general agreement with the 
analysis done in [9] and other data from the literature 
[10, 11, 33]. As already shown in [14] and in general ac- 
cordance with the analysis done in [9], for intermediate 
SPS energies, E^b/A ~ 30-160 GeV, x*/dof is consid- 
erably smaller for the inhomogeneous freeze-out surface 
than for the homogeneous case, which is far outside the 
95.4% confidence interval [35]. At E Lah /A = 20 GeV 
and at RHIC energies, x 2 /dof is similar for the inhomo- 
geneous approach with local strangeness neutrality and 
the homogenous model. However, between 20 and 80 
GeV the x 2 /dof values for the inhomogenous approach 
with local strangeness neutrality are rather large (be- 
tween 2 and 4). In contrast, the inhomogenoues model 
with global strangeness neutrality gives x 2 /dof « 1 for 
E-LaxJA ~ 20 — 160 GeV. It is important to note that this 
result is not due to introducing an additional parameter, 
but just due to allowing for domains of finite strangeness 
with global strangeness neutrality! The calculations us- 
ing global strangeness conservation for RHIC energies are 
under way, but due to the corresponding small baryon 
chemical potentials at these high energies no consider- 
able effect should be expected. Thus, the inhomogenoues 
model allowing for domains of finite net strangeness gives 
a very satisfactory description (x 2 /dof sw 1) of the ex- 
perimental data for particle abundance ratios from lowest 
SPS energies up to highest RHIC energies. However, at 
RHIC the homogeneous approach already gives a good 
description of the data and the inhomogeneous model 
does not provide a statistically significant improvement. 
Thus, the assumption of a nearly homogeneous decou- 
pling surface can not be rejected there. On the other 
hand, the considerable improvement of the description 
of the data for E^^/A ~ 30 — 160 GeV indicates that 
at intermediate and high SPS energies, the experimental 
data favor an inhomogeneous freeze-out surface. For the 
SPS 20 GeV data the situation is not clear: there is cer- 
tainly a reduction of the x 2 /dof in the inhomogeneous 
approach, but also the homogeneous model gives a much 
better value than for higher SPS energies. Here more 
experimental data are necessary to clarify the picture. 
It is worth noting that in general, the improvement due 
to the inhomogenoues decoupling surface is not driven 
by one single species; rather, the inhomogeneous model 
describes nearly all multiplicities better than a homoge- 
neous decoupling surface [36]. 

To illustrate the significance of inhomogeneities differ- 
ently, we show contours of x 2 /dof in the plane of ST, 
SpB in figs.2, 3, and 4. Here, T and ~p B were allowed to 
vary freely such as to minimize x 2 at each point. Fig. 2 
shows that at RHIC energy, x 2 is very flat in both di- 
rections. This shows again that with the present data 
points, a homogeneous freeze-out model appears to be a 
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reasonable approximation at high energies. In contrast, 
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FIG. 2: x 2 /dof contours in the 5T, 5/j,b plane for top RHIC 
energy, (^/sjvjv = 200 GeV). The other two parameters (T, 
Ji B ) are allowed to vary freely. The x 2 /dof minimum is indi- 
cated by the cross. 

fig. 3 shows that x 2 is relatively flat along the S/j,b di- 
rection, while ST is determined more accurately and is 
clearly non-zero. In general we find that in the approach 
with local strangeness neutrality there is little correla- 
tion between ST and SfiB and that about the minimum, 
X 2 is rather flat in S^b direction for all energies. Finally, 
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FIG. 3: Same as fig. 2 for top SPS energy (£ L ab = 158 GeV) 
with local strangeness neutrality. 

fig. 3 shows the contours at SPS 158 for the case of global 
strangeness neutrality. Now, the x 2 determines the S^lb 
more accurately, favoring relatively large finite values. 
For ST, again, values different from zero are strongly fa- 
vored, which, however, turn out to be generally a little 
bit smaller than in the local f s = case. The better 
description of the data and the better accuracy in deter- 
mining the width of the ^^-distribution can be explained 
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FIG. 4: Same as fig. 2 for top SPS energy (£ Lab = 158 GeV) 
with global strangeness neutrality. 



as follows: If vanishing net strangeness is demanded in 
each single domain, in regions with high baryon chemical 
potential the strange chemical potential has to be small 
to guarantee f s = 0. Thus, the possible increased pro- 
duction of strange particles in domains with high baryon 
chemical potential is restricted and results in the shown 
flatness of the x 2 distribution in SfiB direction. In con- 
trast, if the net strangeness vanishes globally, in domains 
of high chemical potential resulting from a large width 
5/ib,& large number of strange particles can be produced. 
Thus, the x 2 should be much more sensitive to the value 

Of <5/i_B. 

As already discussed in section , an inhomoge- 
neous freeze-out surface or finite values for the width- 
parameters result in different mean emission temper- 
atures and chemical potentials for different particle 
species, c.f. eq. 9. These are shown in fig. 5 for the 
case of local strangeness neutrality and in fig. 6 for the 
case of global strangeness neutrality at selected energies 
in the CERN-SPS range. For the cases shown, the in- 
homogeneities determined from the fits to the particle 
abundances are large. Note that the different values for 
these mean emission temperatures and chemical poten- 
tials result from the convolution of the distribution func- 
tion for a given particle species with the Gaussian prob- 
ability distribution determined by the four parameters 
T, TTg, ST, S^b- For both cases, the effect of the inhomo- 
geneities is evident. For example, anti-protons are typ- 
ically emitted from regions with lower baryon-chemical 
potential than protons; also, heavy particles are concen- 
trated in "hot spots" while light pions are distributed 
more evenly throughout the decoupling volume etc. [37]. 
Figures 5 and 6 also show the differences in the result- 
ing mean emission temperatures and chemical potentials, 
depending on whether local or global strangeness neutral- 
ity is adopted: In the case of local strangeness neutrality, 
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FIG. 5: Freeze-out temperatures (T); and chemical potentials 
{(is), of various particle species at E-^^/A = 40, 80, 158 GeV 
for local strangeness neutrality. 
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FIG. 6: Freeze-out temperatures (T); and chemical potentials 
(Ms}i of various particle species at E^b/A — 40, 80, 158 GeV 
for global strangeness neutrality. 



the emission chemical potentials of baryons and the cor- 
responding anti-baryons do differ much less than in the 
case of globally vanishing net-strangeness. For example 
the mean emission baryon chemical potential for the O 
and the f2 are nearly identical for f s = locally, while 
they are widely separated for the global constraint. This 
results from the above discussed effect of the adjustement 
of the strange chemical potential in the local case. There, 
/is is so small (large negative value) in regions with large 



Hb, that the resulting chemical potential of f2 and fl are 
similar. This is not the case anymore if the strange chem- 
ical potential is determined globally and constant for the 
different domains. Then, the mean freeze-out points for 
the different particle species are spread over a much wider 
range in the T— /i^-plane. On the other hand, the spread 
in temperature is somewhat larger in the local case than 
for global strangeness neutrality, resulting from the larger 
best fit values for the width parameters ST. 



SUMMARY AND OUTLOOK 

In summary, we have shown that inhomogeneities on 
the freeze-out hypersurface do not average out but reflect 
in the ( or midrapidity) , single-inclusive abundances 
of various particle species. This is due to the non-linear 
dependence of the hadron densities Pi(T, /i^) on the lo- 
cal temperature and baryon-chemical potential. Conse- 
quently, even the average 'p i probe higher moments of the 
T and [Ib distributions. In [14] we showed that an inho- 
mogenoues freeze-out model with local strangeness con- 
servation strongly improves the description of the data 
at medium and top SPS energies compared to the ho- 
mogeneous freeze-out. Here we showed that inducing 
global strangeness neutrality, results, without adding an 
additional parameter, in a further reduction of the x 2 
at SPS energies. With the resulting xjdoj w 1 for the 
whole range - from lowest SPS to highest RHIC ener- 
gies. Furthermore, while for local strangeness neutral- 
ity we observed a rather flat x 2 /do/ in S^b direction, 
this is determined more accuratrely if strangeness neu- 
trality is ensured only globally. Rather in this approach 
a high statistical significance for a finite width of the 
distributions for temperature and baryon chemical po- 
tential at medium and high SPS energies is observed. In 
addition we showed how in this region the mean emis- 
sion temperature and chemical potential vary for differ- 
ent particle species. Our results also show that there are 
some characteristic differences in the distribution of the 
resulting mean emission values, depending on whether 
strangeness neutrality is fulfilled locally or globally. Es- 
pecially the separation in the mean emission chemical 
potential between baryons and the corresponding anti- 
baryons is strongly influenced by the adopted strangeness 
neutrality condition. 

Inhomogeneities could also affect the coalescence prob- 
abilities of (anti-) nucleons to light (anti-) nuclei, which 
are also sensitive to density perturbations [39]. Other 
signals, such as two-particle correlations [8, 40], could 
also be analyzed in this regard. Future studies should 
shed more light on whether these inhomogeneities can in- 
deed be interpreted as fingerprints of a first-order phase 
transition. Eventually, one would want to establish more 
quantitative relations between the amplitudes of the T, 
fiB inhomogeneities and the properties of the phase tran- 
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sition, e.g. its latent heat and interface tension. 

Data from GSI-FAIR, the low energy program at RHIC 
and and CERN-LHC will provide additional constraints 
for the evolution of chemical freeze-out with energy. 

To improve the quality of the statistical fits, more data 
on hadron multiplicities would be helpful, in particular 
at the lower end of the CERN-SPS energy spectrum and 
at RHIC. This includes estimates of multiplicities of un- 
stable resonances (p, K* , w, A ...) at chemical freeze- 
out [41]. Data from GSI-FAIR and CERN-LHC will pro- 
vide additional constraints for the evolution of chemical 
freeze-out with energy. 
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